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Abstract. We study weakly-compressive high-frequency plasma waves which are 
superposed on a large-amplitude Alfven wave in a multi-fluid plasma consisting of 
protons, electrons, and alpha particles. For these waves, the plasma environment 
is inhomogeneous due to the presence of the low-frequency Alfven wave with large 
amplitude, a situation that may apply to space plasmas such as the solar corona 
and solar wind. The dispersion relation of the plasma waves is determined from a 
linear stability analysis using a new eigenvalue method that is employed to solve the 
set of differential wave equations which describe the propagation of plasma waves 
along the direction of the constant component of the Alfven wave magnetic field. 
This approach allows one to consider also weak compressive effects. In the presence 
of the background Alfven wave, the dispersion branches obtained differ significantly 
from the situation of a uniform plasma. Due to compressibility, acoustic waves are 
excited and couplings between various modes occur, and even an instability of the 
compressive mode. In a kinetic treatment, these plasma waves would be natural 
candidates for Landau-resonant wave-particle interactions, and may thus via their 
damping lead to particle heating. 



1. Introduction 

The solar wind is the classical paradigm of a nonuniform plasma which is structured 
in space and time (Schwenn and Marsch, 1991) and permeated by fluctuations on a 
wide range of scales. It reveals turbulence and carries large-amplitude Alfven waves 
which mostly originate in the solar corona (for reviews of the magnetohydrodynamic 
turbulence see, e.g., Tu and Marsch (1995); Bruno and Carbone (2005) and the 
plasma kinetics see Marsch (2006)). The plasma waves in the solar wind are thus 
riding on a varying background, and throughout their passage from the corona 
the solar wind particles are observed to evolve non-adiabatically, and therefore are 
supposed to undergo continuous heating out to 1 AU (e.g. Kasper et al., 2008) and 
beyond. It is believed since a long time (Coleman, 1968; Belcher and Davis, 1971) 
that wave-particle interactions (Marsch, 2006) with the ubiquitous fluctuations in 
the electromagnetic fields are the main driver of this heating process. Here we study 
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weakly-compressive high-frequency plasma waves which are superposed on a large- 
amplitude Alfven wave in a multi-fluid plasma consisting of protons, electrons, and 
alpha particles like in the solar wind. The term 'high-frequency wave' refers to the 
wavenumber regime around the inverse ion inertial length £p in contrast to the 
non-dispersive low-frequency MHD limit. 

As is well known, there are two important kinetic resonances that can lead to dis- 
sipative heating of a plasma by wave-particle interactions with waves that propagate 
parallel to the background magnetic field. The first is Landau resonance requiring a 
parallel wave electric field, and the second cyclotron resonance (e.g. Akhiezer et al., 
1975; Hollweg and Isenberg, 2002) that couples to the perpendicular wave electric 
field. The condition for Landau resonance is given by kv\\ — w = 0, where k denotes 
the parallel wave number, v\\ the particle velocity in the direction parallel to the 
background magnetic field and lo the wave frequency. This effect can lead to parallel 
heating of the particles (Lehe et al., 2009). In a low-beta plasma such as the solar 
corona, it is difiicult to fulfill this condition if the wave phase speed is close to the 
Alfven speed, because Vth ^ Va (Chandran et al., 2010). 

The cyclotron resonance is connected with the transverse electromagnetic field, 
and the resonance condition is given by fcf || — u; — nClj = with an integer n and 
the particle gyrofrequency ilj (only n = 1 is allowed in the case of parallel propaga- 
tion). This resonance can lead to pitch-angle diffusion which is indeed observed in 
solar wind protons (Heuer and Marsch, 2007). The perpendicular fiuctuations must 
have high wave numbers in the range of the gyroradius to fulfil the resonance con- 
dition. Alfven/ ion- cyclotron waves (A/ICs) are possible candidates for waves that 
can undergo this kind of interaction. However, their origin and evolution are not 
fully understood, even though their existence in the solar wind was recently proven 
(Jian et al., 2009). Indirect evidence for Alfven-cyclotron heating was already re- 
ferred to by several authors beforehand from simulations (Gary et al., 2005) and 
proton in-situ observations (Marsch and Tu, 2001; Kasper et al., 2008). 

The role of weakly compressive waves in the different phenomena of plasma heat- 
ing is currently under wide discussion (Tu and Marsch, 1994; Chandran, 2005; Bale 
et al., 2005; Kellogg et al, 2006; Chandran et al, 2009; Verdini et al., 2010). Es- 
pecially, kinetic Alfven waves (KAWs) have come into the focus of the debate, 
because they are both transverse and compressive. However, several problems also 
arise from this interpretation of solar wind fluctuations, especially at high wave 
numbers (Podesta et al., 2010). 

It has also been known for a long time that compressibility plays a major role 
in the context of the parametric instabilities of large-amplitude waves (Galeev and 
Oraevskii, 1963; Goldstein, 1978; Lashmore-Davies and Stenflo, 1979; Stenflo and 
Shukla, 2007). These instabilities are always connected with compressive compo- 
nents of the daughter-wave products. Our subsequent analysis is based on the pre- 
vious derivations by Marsch and Verscharen (2011), which treat the density fluctu- 
ations in terms of the longitudinal electrostatic field and a ponderomotivc electric 
field. The connection between density fluctuations and ponderomotive forces in 
the context of parametric instabilities was discussed before by Sharma and Shukla 
(1983) for many different wave modes. However, the work of those authors was fo- 
cused on frequencies around the upper-hybrid frequency, which is much higher than 
the frequencies considered here, and thus beyond the scope of this paper. Electro- 
magnetic circularly polarised waves can also be generated from the high-frequency 
side. For example, Murtaza and Shukla (1984) discussed how an upper-hybrid wave 
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can generate such waves in a two-fluid model. For this purpose, electrostatic ef- 
fects and ponderomotive forces had to be included in the model describing the 
high-frequency pump wave. The importance of electrons in the decay of compres- 
sional Alfven waves was discussed more recently by Brodin et al. (2008) in terms 
of the Hall-MHD description. These authors found a new decay channel for oblique 
daughter waves, and discovered that the wave decay products could grow on scales 
around the ion inertial length. They also discussed the role of kinetic Alfven waves 
as decay products and their possible ability to heat the plasma. 

Our present work concentrates on purely parallel wave propagation, but we con- 
sider the role of both electrostatic and electromagnetic components of the wave 
modes in a multi-fluid plasma. The nonlinear coupling of the longitudinal electro- 
static fleld and of the ponderomotive electric fields with the transversal electromag- 
netic wave fields is the main reason for the significant changes we found in the mode 
structure and polarisation. However, the modified wave modes and the possible de- 
cay products are, due to the given geometry, still forced to propagate along the mean 
fleld which is determined by the constant longitudinal magnetic-field component. 
So we do not consider genuine oblique wave propagation. 

The multi-fiuid wave equations are solved here with an eigenvalue and eigenvector 
method that opens a new and unusual way of analysing the pump wave decay and 
dispersion properties of the resulting plasma waves. This approach also provides 
the dispersion and polarisation properties of these waves in a comprehensible and 
direct way. 

As mentioned above, it has been known for a long time that the solar wind is per- 
meated by waves and structures on many diff'erent scales (Tu and Marsch, 1995). 
But also in the solar corona low-frequency waves in the magnetic field were re- 
cently observed by remote-sensing techniques (De Pontieu et al., 2007). Therefore, 
in theory and modelling it is obviously necessary to assume an inhomogeneous back- 
ground magnetic fleld. The natural choice for such a field might be one consisting 
of low- frequency Alfven waves. 

Consequently, the scenario assumed for our present theoretical treatment is the 
following. A low-frequency Alfven-cyclotron wave is assumed to provide the nonuni- 
form background magnetic fleld and corresponding background velocity fleld, ac- 
cording to the wave polarisation relation. In the flank of this wave, a linear disper- 
sion and stability analysis is performed for a three-component plasma consisting of 
protons, electrons, aiid alpha particles, whereby also drifts among these species with 
respect to each other can be included. The situation is sketched in Fig. 1. The ap- 
plied multi-fluid model allows to consider transverse waves with an intrinsic weakly 
compressive component. If transverse dispersion branches at wave numbers close 
to the ion gyroradius are found, then these can be made responsible for possible 
perpendicular ion heating. Their compressive electrostatic components can in turn 
explain parallel heating. Since the observations indicate large perpendicular tem- 
perature anisotropies in the corona (Antonucci et al., 2000; Kohl et al., 2006) and 
in fast solar wind (Marsch, 2006), the assumption of weak compressibility seems to 
be justified empirically, and is also consistent with the measured density fluctuation 
level (Tu and Marsch, 1995). 
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Figure 1. Geometry of the scenario. On a constant magnetic field along the axis of 
propagation fco, a low- frequency Alfven wave (LFAW) propagates and generates the trans- 
verse background fields Bto and Vjto = — C-^to with C, — const. > 0. Weakly-compressive 
high-frequency waves (WCHF) are superposed in the flank of the Alfven wave. 

2. Theoretical approach and numerical treatment 

2.1. The multi-fluid model 

Following the detailed derivations of Marsch and Verscharen (2011), we obtain 
subsequently a set of second-order ordinary coupled differential equations (DEs) 
of the relevant fields for parallel propagation, in which case is a conserved 
quantity, and thus the fields can be made dimensionless as follows: Vjt — Vjt/VA, 
and bt = Bt/B^. Similarly, we define the normalised longitudinal electric fields as 
follows: 

_ cEja _ cEn 



B,,Va' ^" B^Va' ^^-^^ 
Length scales can be all normalised in units of the ion inertial length scale £p = 
c/ujp. In the following, all velocities are normalised to the proton Alfven speed 
Va — Bn/ y/4:TTnpmp and then denoted by lower case letters. This choice is different 
than the normalisation by Marsch and Verscharen (2011) but it is more appropriate 
here, because the wave phase speed V will be used as a free parameter for different 
evaluations of the wave equations. Therefore, it is better to normalise with a fixed 
velocity. The natural choice for such a velocity is the proton Alfven speed. In this 
normalisation system, frequencies are normalised to the proton gyration frequency 
flp, which is equal to Va/^p- In the subsequent paper, we use flj as the normalised 
gyration frequency of species j. We also recall the original definition 

ejn = ^ {vjt X bt) • n. (2.2) 

Using the same normalisation for the transverse electric field, we obtain 

et = ~v{n X bt), (2.3) 

which is fully determined by the solution found for bt, which means it is a depen- 
dent auxiliary field. Completing the required definitions, we quote the differential 
operator 

where we have defined the spatial coordinate along B^ as = which we will 
continue using subsequently. After these preparations we can state the DEs, as 
derived from the equations for parallel wave propagation. Firstly, we obtain for each 
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particle species with respect to its transverse motion a forced harmonic oscillator 
equation reading 

/ e„ + Cjn dt;jt bt , ^ dbt\ ^ 
Vjt = -2 2 — 17- + — +nx — ]flj. (2.5) 

Secondly, one can rewrite the mutually coupled and driven wave equations for the 
transverse magnetic field, 

= - E i f + %^(- X -.t)) , (2.6) 
and the longitudinal electric field, 

^i^e„ = E%- (2.7) 
i J 

All involved parameters such as Vjrv, Cj,\j,ij and rj have non- vanishing mean val- 
ues, as the density nj is always nonzero and its fluctuations throughout are assumed 
to be small. Note that Clj is strictly constant as it depends on the conserved quan- 
tity Bn ■ The other parameters can be calculated by use of the background number 
density fij, as well as by exploiting the conditions for quasi-neutrality, zero lon- 
gitudinal current, and zero centre-of-momcntum velocity (Marsch and Verscharen, 
2011). Compression is accounted for solely by the longitudinal electric field in this 
approximation. This is the major advantage of this multi-fluid system. 

Furthermore, the multi-fluid approach incorporates the natural scales of the 
plasma and allows therefore the treatment of dispersive waves in the high-frequency 
range which is not accessible by MHD considerations. The compressibility can be 
accounted for in a comparably lucid way through the electric field variables Cn and 
ejn- The nonlinear equations (2.5), (2.6), and (2.7) form a closed set of DEs, which 
describe the leading-order compressive eff'ects via en(C)- Before we write down these 
field equations in their components, let us define an appropriate coordinate system. 
For a right-handed orthogonal system, we choose the unit vectors n = 63 = ei x 62, 
ei = 62 X h, and 62 = n x ei . The normalised (dimensionless) transverse magnetic 
field is 

bt = ?'i(e)ei+&2(0e2. (2.8) 
Similarly, the transverse flow velocity of any species j is given by 

Vjt = Vji{£,)ei +Wj2(0e2- (2.9) 

The particle speed along the direction of propagation in the fixed coordinate system 
(not in the co- moving frame) is denoted by Ujn- If there are no differential drift 
motions along the mean magnetic field in the background plasma, i.e., if for all 
j we have Uj^ = 0, then Vj^ = —v, and thus the longitudinal gyration length 
becomes Vj = —v/flj, which by definition is not a positive-definite quantity as 
the gyrofrequency carries the sign of the charge of the species considered. Using 
Eq. (2.8) and Eq. (2.9), the longitudinal electric field associated with species j 
simply reads 

ejn = Vjib2 - Vj2b\. (2.10) 

We recall from Marsch and Verscharen (2011) that e^n can also be expressed as the 
gradient of a potential, which is given by the transverse kinetic energy of species j 
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and then has the normalised form, 



1 d /I 



(2.11) 



In terms of components, we have from a comparison of that equation with Eq. (2.10) 
the relation, 

1 dVji,2 



±b- 



2,1, 



(2.12) 



which we can use to replace the first derivative of the transverse velocity by the 
magnetic field. Written out explicitly, the set of DEs in dimensionless form finally 
reads as follows: 



bi,2 db2,i e„ 



dC^ 



bi,2 ^ 1 

^2 - 2^f2 



711 

— Vj2,l 



-^-^ = T.j2i-nb2-Vj2b,) 



(2.13) 
(2.14) 
(2.15) 



The operators on the left-hand sides of Eqs. (2.13), (2.14), and (2.15) describe 
the dynamics of the uncoupled free fields and have a simple physical interpretation. 
After Fourier transformation (yielding d/d^ ik) the solution of Eq. (2.13) gives 
a helical gyration in space of the transverse velocity about the mean field Bn with 
wave vector k = ±l/rj, the solution of Eq. (2.14) corresponds to diamagnetism, 
i.e. the static penetration of the transverse field into the plasma by the skin depth 
is = c/wp, and Eq. (2.15) gives for t; = the static screening by the Debye length 
Ad, or for finite speed v the electrostatic wave dispersion relation, {kXu)^ + 1 = 0. 
This, for zero drifts, transforms into the electrostatic dispersion relation 



(2.16) 



the zeroes of which which yield the Langmuir and ion-acoustic waves. For finite 
right-hand sides of the above DEs, the fields are coupled (note that the abbreviation 
from Eq. (2.10) has to be included). The incompressible limit (with = e^n = 0) 
then gives the monochromatic (only a single k = k{V) is permitted from the dis- 
persion relation) electromagnetic Alfven-cyclotron wave that has a constant but 
arbitrarily large amplitude. When considering compressibility, the electric and elec- 
tromagnetic waves are linked and interact through the nonlinear rightmost terms 
in the above wave equations. Considering compressibility may require either a per- 
turbative approach or numerical treatment. 



2.2. Perturbative approach and linearisation 

The above equations are a system of second-order wave equations. To simplify it, we 
reduce it to a system of first-order equations. Therefore, we introduce the quantities 

ejj = deii/d^, b'l 2 ^ d&i^2/d^, and w^i 2 ^ duji.2/d^. With this substitution, the 
system corresponds to a set of 18 coupled nonlinear first-order ordinary differential 
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equations in ^ for an electron-proton plasma. For each additional particle species, 
four coupled equations are added. 

In the next step, the system (2.13-2.15) is linearised around a backgromid given 
by the wave amplitude vectors VjtQ and bjto to determine the wave dispersion. A 
suitable choice of the background values is given in Sect. 2.3. 

Nonlinear couplings of low-frequency waves with fluctuations at high frequencies 
cannot be described by a linearised system. The low-frequency wave with wave 
number fcg, however, may be treated as a constant inhomogeneous background if 
the high-frequency waves with wave numbers k fulfill the condition ko <C k. The 
high-frequency waves are then treated in the flank of this low-frequency wave that 
does not change its fields significantly over several periods of the high-frequency 
waves. The particle velocities, the electromagnetic field, and their derivatives are 
combined in a state vector, 

y = iv'pi,Vp2, ■ ■ ■ ,Vpi,Vp2, . . . ,b[,b'2,h,b2,e'^,er,), (2.17) 

and therewith one can write the linearised equation as, 

^5y = A5y, (2.18) 

with the quadratic matrix A = (oj j ) G M{n) with n = 4s + 6 (s is the total number 
of species). The solution of Eq. (2.18) is in general given by 

4s-|-6 

5y=J2 ^iSVie^'^, (2.19) 

i=l 

with the eigenvalues Aj and the corresponding eigenvectors Sy^. Since .4. is a real 
matrix, the complex conjugated eigenvalues and eigenvectors are also solutions once 
a complex eigenvalue or eigenvector is found. The coefficients are arbitrary. 
However, they must be equal for the pairwise complex conjugated eigenvalues to 
construct a real solution. 

An imaginary part of an eigenvalue always indicates a periodic fluctuation. 
Due to the rule of pairwise complex conjugated eigenvalues, they are represented 
by real sine or cosine functions. Real parts of Aj correspond to growth or damping. 
We denote the imaginary part as k and the real part as k. If the eigenvalues are 
pairwise symmetric in the real part, the growth is described by hyperbolic sine or 
cosine functions, which is the symmetric solution for the instabilities growing in 
positive and in negative ^-direction. 

The eigenvalues can be calculated numerically with the QR-method after trans- 
forming A to an upper Hesse matrix (Press et al., 1992). The existence of the 
complex-conjugated eigenvalues means that for wave-like daughter products always 
two solutions exist, one foreward and one backward propagating with the same 
frequencies. 

2.3. Background wave 

An adequate background to calculate the dispersion is the flank of a circularly 
polarised Alfven wave, since these waves are exact eigenmodes of a plasma and 
obey Eqs. (2.13) and (2.14) in the incompressible limit with arbitrary amplitudes. 
The phase speed of this wave is denoted by vo- It is left-hand polarised and has a 
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magnetic field of the structure 

Without loss of generality, we can evaluate this field at the point ^ = 0. For suf- 
ficient low fco-values, this field appears as a constant magnetic background field of 
the magnitude feto = as stated above. Therefore, this background describes 
analogous conditions as used to treat oblique propagation of linear modes. How- 
ever, the wave is additionally associated with a transversal velocity for each species. 
Corresponding to this magnetic field, a background velocity field occurs that is de- 
termined by the polarisation relation of circularly polarised Alfven waves 

as it is derived by Sonnerup and Su (1967) for example. The normal velocity com- 
ponent of the particle species j in the reference frame moving with vo is denoted 
by Ujno- The polarisation relation is also in agreement with the results by Marsch 
and Verscharen (2011). For a sufficient small wavenumber ko, the wave fulfills the 
Alfvenic dispersion relation vq —Va- 

The transformation of the transverse velocity to the co-moving reference frame 
does not change the value for Vjto. Therefore, the polarisation relation provides the 
necessary (constant) value for vjto depending on a given (small) ko and the wave 
amplitude b. 

To evaluate the polarisation relation exactly, the dispersion relation for the cir- 
cularly polarised Alfven waves is used. It is given by 

^o^ + E4ir^^=0 (2-22) 
as the non-relativistic limit from Sonnerup and Su (1967). 



3. Dispersion of high-frequency waves 

In this section we study the wave dispersion relation of high-frequency waves prop- 
agating on the nonuniform background Alfven wave. The most important free pa- 
rameter of the system of equations (2.13-2.15) is the disposable wave phase speed 
V. For a fixed v, the eigenvalue/eigenvector method provides the corresponding 
wavenumber values k{v). The relation co = kv then delivers the corresponding co, 
and thus by scanning throiigh all relevant values of v the hill dispersion relation can 
be determined. The growth rate can be determined similarly by evaluating ^ = kv 
with the spatial growth rate k from the eigenvalue determination. The very high- 
frequency branches, which are dominated by the electron dynamics, are not treated 
in detail, since the ions carry the main momentum, and their wave-induced motion 
is therefore more important for the heating and acceleration processes. The electron 
density and relative velocity with respect to the protons and the other ionic species 
is determined by the requirement of vanishing charge density (quasineutrality) and 
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vanishing constant longitudinal currents of the background, such that 



s 



(3.1) 



S 



(3.2) 



i=i 



The specific- heat ratio is set to its adiabatic value oi jj =5/3 for all species. 
3.1. Homogeneous background 

Before we discuss the case of an inhomogeneous plasma, we apply the system of 
equations and the linearisation to an electron-proton plasma for a homogeneous 
background and with /3j = 0. The plasma beta is defined as the ratio of the 
thermal energy density to the magnetic field energy density for each species. The 
background field is set to -Bn = 5 x 10~^ G and the particle number densities to 
np = Uc = 5cm~^. The ratio B^/up must be parameterised even in the dimen- 
sionless normalisation of the system, in order to fix the normalised plasma fre- 
quency u)j/flp for all electrostatic processes. The result of the calculation is shown 
in Fig. 2a. The two incompressible branches correspond to the transverse particle 
motions. They represent a free motion with the gyration frequency for the protons 
at a; = 1, their normalised gyrofrequency. The branch approaching the gyrofre- 
quency asymptotically corresponds to the Alfven/ion-cyclotron wave. In Fig. 2, the 
cold dispersion relation for Alfven/ion-cyclotron waves 



is additionally plotted as green dotted line. The branch at low fc-values with high 

phase speed v corresponds to the whistler mode shown as red dashed line. The cor- 
rected low- frequency dispersion for the multi- fluid R-mode wave (i.e., the whistler 
wave with ionic effects) is given analytically by 



Both relations show a perfect agreement with the numerical calculations made with 
the eigenvalue method (Stix, 1992). This confirms the validity of our approach. 

In Fig. 2b, the dispersion relation is shown for a plasma consisting of electrons, 
protons, and alpha particles without any relative drifts and without a background 
wave field. The particle species have a beta of 0.01. The proton number density 
is set to rip = 5cm~"^, the alpha particle number density to Ua = 0.04np and 
the electron number density according to Eq. (3.1). The alpha-particle cyclotron 
branch approaches asymptotically the frequency w = 1/2 as expected. This value 
corresponds to the alpha particle gyrofrequency in normalised units. The cyclotron 
branch of the protons is slightly deformed in the low-fc range in comparison to the 
cold plasma dispersion relation. Also the whistler wave mode is slightly shifted. The 
presence of the alpha particles is responsible for this deviation. Further correction 
terms would be needed to represent the dispersion relation in this three- fluid plasma 
analytically. 

Our model also allows to include relative drifts of the particles along the wave 




(3.3) 




(3.4) 
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Figure 2. a) Dispersion relation for an electron- proton plasma in a homogeneous back- 
ground with 13 = 0. The cold dispersion relation for Alfven/ion-cyclotron waves and for 
ion-corrected R-mode waves are additionally plotted, b) The same for a plasma consisting 
of electrons, protons, and alpha particles with /3 = 0.01. 

normal direction. The relative drift speed between protons and alpha particles is 
defined as = Vp^^ — Van- In Fig. 3, the dispersion relation for = 0.2 is shown. The 
alpha-particle branch starting at w = is turned into the so-called beam-mode 
branch (Hne inclined to the left) derived from the resonance condition 

CJ = + kUan, (3.5) 

where Uan — v ~ Va^ denotes again the alpha-particle bulk speed component in 
the wave normal direction in the proton rest frame. For a vanishing drift u^a, this 
mode is flattened back to the horizontal line oj = fla = 1/2. 

3.2. Inhomogeneous background 

Next, our system of coupled wave equations can be applied to an inhomogeneous 
background plasma, which corresponds to realistic solar wind conditions. The plasma 
consists of three species: protons, electrons, and alpha particles. The relative drift 
between protons and alphas is set to Wd = Wpn — fan = 0.2, and the plasma beta to 
= 0.01 for each species. The background wave is assumed to have a normalised 
wave number of k^ = 0.01 and an amplitude of 6 = 0.1. Such a wave has a phase 
speed of almost the local proton Alfven speed. The wave number is small enough 
to allow us to neglect any direct nonlinear couplings between the background wave 
and the high-frequency waves. For the latter, the background wave appears as a 
quasi-constant field with respect to which the system can be linearised. The results 
are shown in Fig. 4. 

The dispersion branches are deformed for several reasons. The proton Alfven/ion- 
cyclotron branch at lower wave numbers is not further deformed as compared to 
Fig. 2b. Some branches turn at some bifurcation positions into completely different 
dispersion branches. These transitions correspond to mode couplings. For example, 
the faster sound wave (alpha-particle mode) couples with the Alfven/ion-cyclotron 
wave of the protons at high wavenumbers, whereas the ion-acoustic wave couples 
with the drift-deformed Alfven/ion-cyclotron wave of the alpha particles. The R- 
mode couples with the ion gyration at w = 1 but, interestingly enough, not with 
the alpha-particle beam mode. It shows a break at w = 1 and then continues above 
that frequency. 

Two additional linear branches occur with merely constant phase speeds of u w 




0.13 and v ~ 0.03. Such linear acoustic modes as these two were also found by 
Mann et al. (1997) for a warm plasma. The wave phase speed of these modes is 
usually determined by 

vphi = — (3.6) 

TOp 

^'Ph2 = (3.7) 

rua 

for a plasma with Ua <^ ?ip- The first is the so-called ion-acoustic speed, the second 
is the sound speed of the alpha-particle component. In our case, these velocities 
are given by ?;phi = 0.13 and vph2 = 0.23. The first one corresponds perfectly 
with the fast (steeper line) sound wave that we have found. The drift velocity Ud 
of the alpha particles is the reason for the deviation of the second linear mode. In 
the non-drifting reference case, the two velocities match (uph2 — Vd = 0.03). If the 
numerical dispersion code is applied to a plasma without drift, the second sound 
wave branch is directly found at v — wph2 (not shown here). 

The two sound wave modes do not appear in the homogeneous plasma. Remem- 
ber, that the phase speed v is a free parameter in our calculation, and a wavenumber 
is given by the eigenvalues obtained for each v. If a wave is not dispersive, and there- 
fore has a constant phase speed, it is not possible to calculate the full set of possible 
/c- values that belong to this v. The sound waves exist already in the homogenous 
case; however, they become visible only after becoming dispersive due to the mode 
coupling. This type of wave becomes only dispersive for non-zero beta and in the 
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Figure 4. Dispersion relation for a plasma consisting of protons, electrons, and drifting 
alpha particles. The frequency uj is shown in dependence on the wave number k in dimen- 
sionless units. The derived A/IC branch deviates from the cold dispersion branch at low 
and high wave numbers. Mode coupling between the various branches occurs as discussed 
in the text. 

presence of the background wave. Also above the corresponding gyrofrequencies 
CO = 1 and w — 1/2, the branches continue with their constant phase speed. 

The waves on all branches have a compressive component, owing to the non- 
vanishing electric fields en and Cjn- The real part of the eigenvalues Xi is always 
zero, i.e., none of the compressive modes is unstable in this case. 

Note that it is the amplitude of the background wave which mainly determines 
the strength of the mode coupling and, therefore, the position and shape of the 
deformation of the ion-cyclotron branch. For higher amplitudes, the branch turns 
earlier away from the A/IC dispersion branch and, hence, the phase speed of the 
linear mode increases already at lower fc-values. The plasma beta determines the 
phase speeds of the linear branches. Higher betas lead to higher phase speeds of 
these modes. The phase speeds can be adjusted relative to each other by choosing 
different betas for the individual species. However, the overall topology of the dis- 
persion branches is not significantly changed by assigning different beta values to 
the different species. 

However, if beta is chosen to be greater than one, the situation changes tremen- 
dously. The sound waves become unstable. In Fig. 5a, the dispersion is shown for 
a situation with /3j — 1.2. The growth rate of the modes is shown in Fig. 5b. 

Mode coupling between the sound wave and the whistler wave mode leads to 
an instability that occurs at /c ~ 0.8... 1.6. Interestingly enough, this range is 
the ion gyroradius scale for moderate plasma betas close to one. This growing 
instabihty can therefore provide wave energy to the dissipative regime by means 
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Figure 5. a) Dispersion relation for a plasma with /? = 1.2. b) Imaginary part of tlie 

frequency depending on k. 



of Landau resonance, since the particle distribution has a significant number of 
particles at the corresponding wave speed of about v « 1.36 for a plasma with 
(3j > 1. This instability vanishes if the background wave is absent, i.e., in the 
homogeneous case. For higher beta values, the position of the instability is shifted 
to higher wavenumbers. 



4. Conclusions 

As expected, the inhomogeneity of the background wave field leads to a deforma- 
tion of the standard normal modes in the plasma. All of the discovered waves have 
a compressive component. Compressibility and inhomogeneity of the background 
plasma are the reasons for the new effects in the dispersion properties. The initially 
non-compressive Alfven/ion-cyclotron waves become slightly compressive due to the 
inhomogeneous background, and are therefore able to interact with the protons via 
their electrostatic field components, in addition to the cyclotron resonance enabled 
by their transverse components. The quantitative details of these wave-particle in- 
teractions cannot be treated in the fluid description, since they are purely kinetic 
processes. Therefore they require a kinetic Vlasov treatment at least in the frame 
of quasi-linear theory and are, therefore, beyond the scope of this work. The ki- 
netic refinement of the dispersion analysis of circularly polarised waves of the same 
type was described by Stenflo (1976), who also showed how relativistic effects and 
compressibility can modify the dispersion of large-amplitude waves in multi-fluid 
theory. 

The non-constant (due to the presence of the pump wave) background leads to 
nonlinear mode couplings between some of the linear wave modes. It also leads to 
the excitation of initially non-dispersive modes, such as the ion-acoustic wave or 
the alpha-particle sound wave, yet now with fc-dependent wave speeds. 

The acoustic modes can grow at wavenumbers around the ion gyration scale 
under certain conditions and for adequate parameters. Thus, they are good candi- 
dates for a longitudinal electrostatic wave field with which particles can undergo 
Landau-resonant wave-particle interactions. Kinetic Alfven waves (KAWs) have re- 
cently been discussed as another possible reason for this kind of interaction. KAWs 
owe their compressibility to the oblique geometry of their propagation with re- 
spect to a constant background field. In this study, we can show that also waves 
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with purely parallel propagation can grow nonlinearly at the corresponding res- 
onant wavenumbers, provided that the background has a non-trivial — but quite 
reasonable — magnetic field configuration. These two mechanisms are quite differ- 
ent in nature and should be further investigated and compared with each other. 

The instability of ion-acoustic modes coupled to transverse modes were also dis- 
cussed in studies of the parametric decay of largo-amplitude; pump waves and evalu- 
ated by numerical simulations (Araneda et al., 2007; Valentini and Veltri, 2009). In 
these studies, also instabilities of longitudinal waves at high wavenumbers are found, 
however at different background parameters. These waves are usually interpreted 
as results of nonlinear wave-wave interactions. Our linearised wave equations, how- 
ever, can explain a similar growth of daughter waves at higher frequencies than the 
initial pump wave frequency as a consequence of compressibility and a non-uniform 
background. This mechanism can thus be understood as a new spectral transfer 
process of plasma fluctuations. It is very similar to the decay instability, which is 
found with the characteristic fc 3> /cq. 

Our eigenvector analysis still keeps the freedom to choose the amplitudes in 
Eq. (2.19). This means that a wave can in general only occur if its amplitude is finite. 
It is beyond the scope of this work to investigate how the discovered wave branches 
can be excited in a real plasma. The linear dispersion analysis can only show possible 
normal modes. The discovered instability leads to growth in the initial phase only 
until nonlinear couplings and perhaps saturation occur. Yet unstable modes can 
grow from the thermal noise (that is constrained as a finite eigenvector of the system 
of equations) with a certain finite amplitude. In the presence of a large-amplitude 
wave, this thermal noise can lead to growth according to the calculated growth rate 
7 at the gyroradius scale range for /3 > 1. Maybe, full nonlinear calculations can 
investigate the further evolution and possible nonlinear excitations of such modes. 
The linear approximation appears to reflect the basic situation well (Lehe et al., 
2009). 

In the future, also other background conditions should be inquired. The above 
chosen background is one of the simplest inhomogeneous conditions. It is important 
to note that in our approach the background has to change slowly in dependence 
upon the position ^. Otherwise, the use of a fixed phase speed v is not possible 
anymore, and then the present approximation needs to be changed. The original 
set of coupled wave equations remains valid, yet another mathematical treatment 
is required to cope with their nonlinearity. 
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